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Osculatory Interpolation in the Complex Plane 

Herbert E. Salzer' 

Tables of coefficients to facilitate osculatory n-point interpolation (n = 2(l)7) in the complex 

plane are given. 

The writer has described in a previous article [1] 2 a new method for oscillatory interpolation 
that can be applied also in the complex plane when one has an analytic function f(z), z=x-\-iy, 
which is known for arguments z k at equal intervals along any straight line in the 2-plane. 
However, for the tabular arguments z k =x k -\-iy k in the form of a Cartesian grid of length h, 
greater accuracy may be had by basing the Hermite osculatory interpolation formula upon 
values of f(z k ) and /'(£*) at points z k , which are closer together and not necessarily lying upon 
a straight line (except, of course, in the 2-point case). For a detailed discussion of the Hermite 
osculatory interpolation formula, the reader is referred to the previous article [1]. This 
present paper merely aims to supplement that article which was intended primarily for real 
functions, by giving here the corresponding auxiliary quantities, a k and b k) which are suited 
better for complex interpolation. Whereas the quantities, a k and b k , for real interpolation 
were tabulated up to the 11-point case (21st-degree accuracy), this present tabulation for 
complex interpolation does not go beyond the 7 -point case (13th-degree accuracy), which is 
more than adequate for most of the practical problems that would arise. 

To interpolate for f(z), where z=z +Ph and P^p-^iq is now complex, we choose the 
z k =Zo-\-kh for k equal to certain small complex integers and such that z Q will always be the 
lower left corner point of the configuration of points, z and z k . Also, in almost every example, 
!^!<1> |<?!<1> so that the point z is usually not outside the square whose corners are z , z u z u 
and Zi+i. For each of the w-point formulas, n = 2(l)7, the points z k are chosen to lie in the 
following configurations : 



Two 


-point 


Three 


-point 


Four-point 




Five-point 








•z { 




•Zf 


'Zi +i 


% Zi 


'Zi + i 






'Z 


'Z 






'Z 




'Z 




■3d 


•Zi 


•Zo 


'Zi 


•Zo 


•Zl 


% 


*Zi 


•z 2 






Six- 


point 




Seven-point 







'Z 2 i 



•z 2i 



'Zi 



'Zl+i 



'Zi+i 



'Z2+i 



'Z 



'Zi 



•Z 



•0o 



J Present address: Diamond Ordnance Fuze Laboratory, Department of Defense. 
2 Figures in brackets indicate the literature references at the end of this paper. 
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Hermite's n-point osculatory interpolation formula for j{z) = /, in terms of f(z k ) =f k and 
£'(z k )=j' hi is expressible in the concise form 



where 
and 



k k 

a k =aJ(P-k) 2 +h/(P-k) 
fama,J(P-k), 



CD 

(2) 
(3) 



the summation being taken for Jc corresponding to the n-points z k of the configuration. The 
dependence of a k , fc , a k , and b k upon n is understood. The auxiliary quantities, a k and b k , 
are tabulated below for n=2(l)7 for each of the above n -point configurations of points z k . 
The actual definitions of a k and b k are as follows: 



a k =C(n) 



aud 



'{f M }' 



b k = -2C(n) I Jp U'(P-j) y J I W(k-j) X, 



(4) 



(5) 



where the product is taken for j corresponding to the n— 1 points Zj9^z k of the configuration, 
and where for each n, C(n) is chosen as the least integer which allows the quantities, a k and b k , 
to appear as complex integers. 



Table of a t and b{ 
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Formulas (1), (2), and (3) for complex osculatory interpolation, which utilize the Cartesian 
grid, will be especially useful in connection with such tables for complex arguments as (1) 
logarithm of the gamma function, log T(z), together with its derivative the psi function, ^(z). 
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(2) Bessel functions of the first or second kind giving Jo(z) and J 1 (z) = —Jq(z) j Y (z) and 

Y 1 (z) = — Y' (z), or linear combinations of (hem, (3) probability Integra e~ u * du with its 

integrand, (4) miscellaneous tables of integrals of the more elementary functions where the 
first derivative or the integrand, although not tabulated, is easy to calculate, namely, the func- 

r* oo 

tion I (e~ u lv)dn and (5) tables of solutions of important linear differential equations, together 

with their first derivatives [2]. In all such tables, and in many others, the user will find these 
complex oscillatory interpolation formulas to be particularly convenient. They are especially 
suitable in those cases where the grid length h is too large for sufficiently accurate complex 
interpolation, using either tables of complex interpolation coefficients [3, 4] or formulas cor- 
responding to (1), (2), and (3) above [5], where the tables or formulas involve the use of only 
the functional values//. 

Illustration. — To demonstrate the use of these tables in formulas (1) to (3), consider 
an example from [2], where the modified Hankel functions hj(z)=(l z^H { ^(^ 0*), j=l, 2, and 
their first derivatives are tabulated over a Cartesian grid of length h=0A in the complex plane. 
Suppose that it is required to find h x (1.24579 316 + 0.96155 803i), using the four-point 
oscillatory interpolation formula and the tabulated values of /^, h[ at z =1.2 J r0.9i, 2i=1.3 + 0.9i, 
Zi=1.2 + l.0i, 2 I+1 =1. 3 + 1.0? [2, pp. 21, 23]. Employing (2), (3) and a u b t from these tables, 
since P=0.45793 16 + 0.61558 03i, one finds that 

ao=i/(0.45793 16 + 0.61558 03?') 2 +(3 + 3?)/(0.45793 16 + 0.61558 030, 
ai =-i/(— 0.54206 84 + 0.61558 03i) 2 +(-3 + 3i)/(-0. 54206 84 + 0.61558 03?), 
a t =-i/ (0.45793 16-0.38441 97?) 2 + (3-3^/(0.45793 16-0.38441 97*), 
a H .,=?/(-0. 54206 84-0.38441 97i) 2 + (-3-3i)/(-0. 54206 84-0.38441 97i), 
O =?7 (0.45793 16 + 0.61558 03i), fr = -?/(-0. 54206 84 + 0.61558 03i), 
Pi=-i/ (0.45793 16-0.38441 97?'), ft +*=?'/( -0.54206 84-0.38441 97i), 

from which one obtains. 

«o=7.09824 5538-1.29187 9023?', j3 = 1.04576 6212 + 0.77794 7889i, 

^=6.63654 0823 + 0.51581 3695?, ft= -0.91498 3814 + 0.80571 7487?, 

a,=9.82415 3465-1.10147 8051?, ft= 1.07536 0453-1.28099 9733?, 

a 1+ ,= 8.43081 1380+1.81986 9385?', ft +< = -0.87048 2793-1.22746 3667i, 

Then for (1) we employ the following tabulated values of hi and h[: 

/ =0.19018 585-0.19313 840?, /o=0.14952 614+0.33229 46K, 

fi = 0. 20311 754-0.15966 418i, /J = 0.10901 235+0.33659 23H, 

/,=0.15891 089-0.17847 850?;, f' = Q. 14376 041+0.29375 232?, 

f l+i =0. 17143 575-0.14888 108?:, /i +t . =0.10669 347+0.29764 346?', 

and also 

2>*=31.98975 121-0.05767 399i, 
k 

so that we have 

^(1.24579 316+0.96155 803t)=/~[a /o+«i/i +«^ 

+a 1+ ,/ 1+ ,+o.i{/3o/;+ft/;+ft/;+^ + ,/; + ,}]/s^ 
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whose terms are 

[(7.09824 5538-1.29187 9023i) (0.19018 585-0.19313 840i) 

+ (6.63654 0823 + 0.51581 3695i) (0.20311 754-0.15966 418?) 
+ (9.82415 3465-1.10147 805 K) (0.16891 089-0.17847 850?') 
+ (8.43081 1380+1.81986 9385i) (0.17143 575-0.14888 lOSi) 
+ 0.1{ (1.04576 6212 + 0.77794 7889?) (0.14952 614 + 0.33229 46 K) 
+ (-0.91498 3814 + 0.80571 7487^) (0.10901 235 + 0.33659 231?;) 
+ (1.07536 0453-1.28099 9733^) (0.14376 041+0.29375 232?) 
+ (-0.87048 2793-1.22746 3667$ (0.10669 347 + 0.29764 346?;)}] 
-*- (31.98975 121-0.05767 399i), 

or after multiplying, 

[(1.10047 4414-1.61664 0896?) + (1.43035 4817-0.95484 7040?) 

+ (1.36457 4820-1.92843 7032?;) + (1.71628 6592-0.94319 763 U) 
+ 0.1 {(-0.10213 8505+0.46382 602H) + (-0.37094 2846-0.22014 3359i) 
+ (0.53089 0903 + 0.13173 258H) + (0.27247 1703-0.39005 5868?*)}] 
-H31.98975 121-0.05767 399i) 
-(5.64471 8767-5.44458 6662i)/(31.98975 121-0.05767 399?;), 

or finally, 

0.17676 025-0.16987 916?;, 

which is correct to within a unit in the last place, as was seen from the independent calculation 
of ^(1.24579 316 + 0.96155 803?) from the Taylor series around z Q = 1.2 + 0.9?, 
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